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ABSTRACT 

This four- volume final report is concerned with Phase II of 
the DOE/JPt project "RF Model of the Distribution System As a 
Communication Channel." An earlier Phase I effort was concerned with 
the design, implementation, and verification of a computerized model 
for predicting the steady-state sinusoidal response of radial (tree) 
configured distribution feeders. That work demonstrated the feasi- 
bility and validity based on verification measurements made on a 
limited size portion of an actual live feeder. The Phase II effort 
Is concerned with 1) extending the verification based on a greater 
variety of situations and network size, 2) extending the model capa- 
bilities for reverse direction propagation, 3) investigating para- 
meter sensitivities, 4) improving transformer models, and 5) investi- 
gating procedures/fixes for ameliorating propagation "trouble spots." 


PREFACE 


This volume contains a collection of miscellaneous appendices, most of which are essen- 
tially verbatum copies of selected internal memos on subjects of general interest. The 
subject/titles are listed on the following pages. 
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APPENDIX 1 


TUTORIAL REVIEW OF “PERTURBATION” THEORY AS USED 
IN DIFNAP PROGRAM NTWKANSj 

R.C. Rustay 


Late in the Phase I work, J.T, Gajjar conceived an exact perturbation theory which enables 
efficient calculation of the effects of single load variation (say, for example, short circuit fault 
to ground effect) and propagation due to any single source not located at the original root sec- 
tion source. This work was reported in the following reports by J.T. Gajjar. 

[1] Matrix Analysis of Effects of Load Variation on Signal Propagation on Multicon- 
ductor Lines, 7/14/77. 

[2] Reciprocity in Multiconductor Systems With Application To Effects of Load Per- 
turbations on Multiconductor Lines. 

This Appendix, benefitting from the insight gained from these reports, presents an alter- 
native approach to the understanding of this theory. 

Further relevant comments are contained in Appendix 2 concerning using reciprocity and 
precomputed driving point admittance to compute reverse direction voltage transfer ratio ma- 
trices. 

Circuit Analysis Background 

The following is a summary of pertinent circuit analysis which will be of use in explaining 
J. Gajjar’s Perturbation Theory. Consider first a two port network with equal number** of 
terminals in each port and represented by ABCD parameters, i.e., 



E lt E 2 , I\, / 2 each n-vectors 
Sketch 1 


** A requirement for the use of ABCD representation with invertibility capability. If the number of terminals are 
not equal, then “one way” representation may be possible if “controllability” conditions are satisfied. Even 
with equal number of terminals, “controllability” conditions may not be satisfied. In applications of interest 
bidirectional controllability is assumed. 
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The circuit equations are*; 

* Ei - AE 2 + BI 2 A-l 

for Sketch 1 

I x - CE 2 + DI 2 A- 2 

from which (and using other fundamental circuit properties) can obtain 

E 2 ~ D T E X - B T I X A- 3 

for Sketch 1 

/ 2 = - A T l x A-4 

First note the reciprocity of transimpedance and transadmittance, 

The transadmittance is concerned with the short circuit output current resulting from an 
input voltage. From A-l with E 2 — 0 

E : - 5/ 2 E 2 « 0 A-5 

so that ^~ 1 r is the transadmittance from port 1 to port 2. Similarly from A-3 with E x — 0 

E 2 =*-B t I x E { - 0 A-6 

so that (taking into account the positive sign/direction of /j defined in the Sketch 1) i?“ lr 
represents the transadmittance from port 2 to port 1, and it is again seen that these two trans- 
admittances are related simply by their transpose, i.e,, reciprocity. 

The transimpedance is concerned with the open circuit output voltage resulting from an 
input current. From A-2 with I 2 — 0 

1 1 - CE 2 / 2 = 0 A-7 

so that C -1 is the transimpedance from port 1 to port 2. Similarly from A-4 with I x — 0 

/ 2 - - C T E X /j - 0 A-8 

so that (taking into account the positive sign/direction of I 2 defined in the Sketch 1) C~ lT 
represents the transimpedance from port 2 to port 1, and it is seen that these two transim- 
pedances are related simply by their transpose, i.e., reciprocity. 

It is of interest to note that the transimpedance which can be derived from equations A-l 
and A-2 with I x — 0, i.e., from port 2 to port 1, can easily be shown to be 

Ei- [ 5 - AC~ x l\l 2 A-9 

which implies that 

C~ lT = [ AC~ X D - b\ A-10 

This can be proved utilizing, appropriately, various identities relating the A, B, C, D matrices. 
This exercise illustrates the “redundancy” of the ABCD parameters and that results derived 
using them should always be carefully examined for possible simplifications. 

Next consider the following circuit. 


* The ABCD matrices are related by A T D—C T B =■* /from which many equivalent relationships can be derived. 
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Sketch 2 


Using equations A-l and A-2 the following relationships are easily derived 


Yjy — ( C+DY 2 )U+BY 2 H 

for Sketch 2 

Y out - (D+Y^y^C+Y^A) 


A-l 1 
A-12 


where Y IN is the driving point admittance looking into port 1 with (only) Y 2 connected to 
port 2. Similarly Y 0UT is the driving point admittance looking (back) into port 2 with (only) 
connected to port 1. 


Next consider the equivalent circuit corresponding to Sketch 3, i.e., 



Sketch 3 


It is easy to show using equations A-l and A-2 that the A e , B e , C e , D e matrices for the 
equivalent circuit (bounded by dashed lines) are 

A e ~A+BY 2 A-13 

B e ~ B A-14 

c e — r,A + c + ( r,B+D ) y 2 for Sketch 3 a-is 

D e - y t B+ D A-16 

so that 


E\ — A e E 2 + B e I 2e 

he “ Q-^2 Dehe 


for Sketch 3 


A-17 
A- 18 


and their inverse from A-3 and A-4 










Sketch 4 

Using equations A-9 and A- 10 the following equation is easily obtained 

E 2 = \a+BY 2 \~ 1 \ Y x HC^DY 2 )U J cBY i y x \' X I^ A-21 

voltage transfer Y/^ 

or after some simplification (or directly from A-18) for Sketch 4 

E 2 — \y x U+BY 2 ) + C+DY 2 y l I { A-22 

What has been obtained is the forward (port 1 to port 2) transimpedance for the total net- 
work denoted by the dashed iine boundary. Equation A-22 will be used to develop the fun- 
damental perturbation equations. Note that the result A-21 can be directly obtained using 
A- 5 and A-18, i.e., C~ l . For port 2 to port 1 transimpedance use A-20 

*1 = - C'T-'lu 

as expected. 

Finally consider the following circuit with a drive / 2 
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Using reciprocity stated earlier 

- Y X (A+BY 2 ) + C+DY 2 ]~ lT I 2 A-23 

V * *\ 

It is of interest for subsequent analysis to develop a rather obvious expression for £ 2 * Using 
A-12 obtain 

E 2 - - \ Y 2 +(D+ r!S)->(C+ Y\A)\~'i 2 A-24 

v "'— V — ■ — ■■■■ ■ ' 

r OUT 

which can be expanded to yield 

E 2 = ~ (D+Y x B)Y 2 + (C+Y x A)^\d^ Y X B) I 2 A -25 

This same result can be obtained from A- 18 with I Xe - 0. 

Development of Exact Analysis for Single Load Perturbation 

The procedure for exactly computing changes in voltages due to a change in a (single) 
load was originally devisbd by Prof. J. Gajjar and reported in references [1] and [23. The 
analytical development to follow arrives at the same procedure via a somewhat different logi- 
cal route, which however, has depended on the insight obtained from Gajjar’s approach. 

Frequently we are interested in the following circuit with a current drive I x and zero out- 
put current (i.e., all loading is represented by Y 2 ), 


0 


Sketch 6 

The current I x and voltage E 2 is given by (see A-22) 

I { - [Y X {A+BY 2 ) + C + DY 2 ] E 2 for I 2 = 0 D-l 

Also note the voltage transfer relationship (see A- 17). 

E x » [A + BY 2 ) E 2 for \ 2 = 0 D-2 

Now suppose that the load Y 2 changes from Y 2 to r 2 + A Y 2 while I x and other circuit 
parameters remain unchanged. Then E 2 will change to E 2 +EE 2 and E x will change to 
E x +hE x . By substituting 

E x E x + hE x D-3 

E 2 ~* E 2 + A E 2 D-4 

Y 2 -+ Y 2 + EY 2 D-5 
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into equation D-l and D-2 it will be shown that A and A£ 2 resulting from the change A F 2 
can be evaluated computationally by driving the output side of the circuit in Sketch 6 by a 
current -A F 2 (£ 2 d-A£ 2 )* i.e„ 



* Sketch 7 

and that A£ 2 can be gotten from Af 2 using various circuit parameters associated with the 
original unperturbed networks, 

As a first step make the substitutions D-3, D-4, D-5 into the transimpedance relationship 
of D-l and then subtracting D-l obtain the following equivalent forms, each of which are 


exact 

A£ 2 - —[(£>+ F t £) JVHC+ F^hKZM- F,£)(£ 2 +A£ 2 ) D-6 

A£ 2 « -[ F 2 +(jD+ F 2 (£ 2 +A£ 2 ) D-7 

h)ir 

A £, - - [ lVfA F 2 + (£+ F, £)“ l ( C+ 1 % 4 V HA F 2 £ 2 D-8 

>oir 

Next make similar substitutions into the voltage transfer relationship of D-2 and then sub- 
tracting D-2 obtain the following exact form 

A E x - 5A F 2 (£ 2 +A£ 2 ) + U+£F 2 )A£ 2 D-9 

Substituting D-7 into D-9 obtain the exact form 

A£j — {£ — C4+£F 2 HF 2 + (D+ F^hKC+F^H) A F 2 (£ 2 +A£ 2 ) D-10 


Now make some observations. First compare Sketch 7 with Sketch 5 and correspondingly 
equations D-6, D-7 with A-25, A-24 and conclude that If a current source A F 2 (£ 2 +A£ 2 ) is 
computationally applied to the unperturbed circuit, in the manner shown in the sketch, then 
the voltage at port 2 will be A £ 2 . Similarly using equations A-17, A-18 (with 7je=0) with 
equations A-13, A-14, A-15, A-16, the result of D-10 can be obtained with 

/ 2 e— A F 2 (£ 2 +A£ 2 ) and the voltage at port 1 being A£|. Hence, if this equivalent current 
A F 2 (£ 2 +A£ 2 ) can be evaluated, then the voltages A£ { and A £ 2 can be evaluated by compu- 
tationally driving the unperturbed circuit as shown in Sketch 7, 

Also by noting that the form of equation D-10 is exactly that of A-9, the reciprocity rela- 
tionship of equation A-10 can be used to obtain the relationship (see A-23) 

A£| - - l )\ (A + B Y 2 )+ C+ D F 2 ]- [ r A F 2 (£ 2 +A£ 2 ) D-ll 

as a replacement for D-10 


1-6 





ORIGINAL PAGE ££> 
OF POOR QUALITY 


Finally note that the equivalent drive current A Y 2 (E 2 +k E 2 ) can be evaluated if A£ 2 can 
be computed, which can be done using equation D-8 which involves known and available pa- 
rameters, 


Thus the computational procedure for evaluating the change A£i caused by a change in 
A Y 2 is as follows: 

a) Compute A E 2 using D-8, noting that Y our , Y 2 , Y x are available from the previously 
computed unperturbed network for any two locations comprising port 1 and port 2 

b) Compute the “drive” current A Y 2 (E 2 +kE 2 ) noting again that E 2 is also available at 
every port of interest in the original unperturbed network analysis. 

c) Determine the forward transimpedance from port 1 to port 2 using y, N and the vari- 
ous voltage transfer ratios also available from the original unperturbed network analysis 

d) Transpose the transimpedance found in step c) above and use it in equation D-ll to 
find the desired voltage change AiTi at the input port. 

Specifically then to carry out the exact perturbation analysis for a single load variation all 
that remains is to establish the forward transimpedance from port 1 to port 2. Consider the 
following sketch which snows a portion of a total (tree) network 



Sketch 8 


where it is of interest to determine the change in voltage of Ej (and E k ) due to a change in 
Y k , a load connected to port 2 of section J. As a result of analyzing the unperturbed network, 
the following quantities (vectors and matrices) are available. 

YquTj ■ equivalent Norton’s source admittance looking back into the output port 2 of 
section j— 1 

I N - equivalent Norton’s source current out of section j— 1 
YiNj - admittance looking into port 1 of section j 
E k - original voltage at port 2 of section j 

Y 0 [/T k ' equivalent Norton source admittance looking back into the output port 2 of 
section j 

Y lN k - admittance looking into port 1 of section k=*j + 1 

Yj - external load (not including Y JNj ) connected to port 2 of section j— 1 

Y k - external load (not including connected to port 2 of section j 


1-7 






ORIGINAL FA(*i & ' ! -> 
OF POOR P^ALIfY 


• voltage transfer ratio from port 1 to port 2 of section j. Note this voltage 
transfer ratio may itself be the product of many sub-sections in cascade which 
form section j\ 

Then to conduct the perturbation analysis associated with a change A Y k in Y k , form or com- 


pute 

D “ 12 

^ 2 “ Y k + Y lNk D-13 

E 2 - E k D 14 

Ar 2 ~Ar A . D-15 

- 7^1 Y out + Yj]~ 1 forward transimpedance across section j D-16 

A£ 2 - - [ K 2 +A Y 2 4- Y 0UTi )-i& Y 2 E 2 D-17 

A £, - - C“ lr A r 2 (£ 2 +A£ 2) ~ [r 0(/r + l}]- 1 T/A Y 2 (E 2 +bE 2 ) D-18 

Note that either or both Yj and Y k may be zero (null matrices), 

Development of Exact Analysis for Reverse Direction Propagation 


This procedure is concerned with the computation of voltage propagation in the reverse 
direction, i.e,, considering Sketch 9 


E 1 

A B 

C D 

E 2 



y in^ 

► 

Y 

^ — out 




Sketch 9 


with any relavent load and termination admittances absorbed into the equivalent A, B, C, D 
as shown in Sketch 3 and equations A-13 thru A-20, what is E { given £ 2 ? Letting, as before, 
y 0 ut to be the admittance looking backward into port 2, and using the transimpedance 
reciprocity relationship, then 


reverse direction transimpedance 

' E\ *■ C- [t Y 0UT E 2 E-l 

This is easy to show analytically, i.e., usinft A-l, A-2 with I x AO obtain 

E x - [A-BD~ l C]E 2 - [AC~ l D~B] D~ l CE 2 E-2 

From A-10 note AC~ X D~B**C~ { T the reverse direction transimpedance and for the 
“unloaded” circuit of Sketch 9 ZHC is Y 0UT (see A-12 with Ti=-0), Thus all that is required 
is to develop the transimpedance as explained previously for the load perturbation case. All 
necessary circuit parameters such as Y 0UT , Ty etc are known and available from the original 
unperturbed circuit analysis. 
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The previous form suggests an equivalent Norton source current effect. Consider the fol- 
lowing circuit 



1 


where Y 0 and I 0 represent a Norton’s source. With /j—0 it is easy to show that 

Ei = l (A+BY o )(C+Dr o )~ l D-3]/ 0 E-3 

Now using A-13, A-14, A-15 and A-10 note that the matrix [(A+ BY 0 )(C+ DY 0 )~ X D- B) is 
the reverse direction transimpedance. (This can be observed easily by assuming Y 0 to be 
imbedded in A, B, C, D and letting ■ 0 yielding AC~ X D—B as before. Hence finding the 
equivalent Norton’s source may be more convenient than specifying the voltage E 2 . 

A final note of interest is to expand equation E-l into terms involving external behavior, 
i-e. Y in , Y out , and the voltage transfer matrix A~ l , i.e. 

E 2 — A~ l E { E-4 

Again emphasizing that all load and termination admittances are absorbed into the ABCD net- 
work of sketch 9 then 


y, N - ca-' 

E-5 

so that 


C" ir » y in ~ x A~ xT 

E-6 

Substituting into E-l 


E\ Y [N ~ [ A~ [T Y 0 ut E 2 

E-7 

where A~ l is recognized to be the forward voltage transfer ratio matrix, Y IN , 

Y out and the 


“overall forward voltage transfer matrix” are all available (i.e, computed and saved) and pro- 
vides a compact procedure. Equation E-7 has been derived using ABCD network representa- 
tion which implies A to be square (equal number of terminals at input and output ports). 
However, if 


then 


E 2 “ RE\ 


E-8 


*i- Y in ~ x R t Y out E 2 E-9 

is perfectly general whether or not input and output ports have equal number of terminals. 
This can be shown using general admittance matrix representation of the two port network or 
more directly utilizing the reciprocity of transimpedances, i.e., (with suitable sign conven- 
tions) 

E 2 !aB ^12 Y in E\ =“ > R = Z n Y in =* > Z 21 = Z 12 r * Yj N ~ x R T E-10 
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£| °* Z 2l Yqvt ^2 am Yin 1 R r Yquj E 2 E-ll 

Note on Cumulative Voltage Transfer Matrices 

The preceding discussions have shown how reciprocity properties makes possible the use 
of a cumulative forward transfer function matrix (R) to compute reverse path propagation to 
the original source. Also implied in this approach is the possibility of computing point to 
point propagation. However, in using reciprocity and cumu’ative voltage transfer matrices for 
this purpose, a subtle computational situation arises as pointed out by J.T. Gagjar. This prob- 
lem can be illustrated by the following sketch. 



Sketch 11 


The voltage transfer ratio matrix T t j relates input and output by 


Ej = T U E, F-l 

Suppose that this network has been analyzed in the conventional outbound direction so 
that the voltage transfer ratios T\ 2 r T 13 , and 7 ^ 4 * are available (along with the various Y lN 
and Y 0 ut admittances), and it is desired to propagate from point 3 to point 4, or from point 3 
to point 2. In either case it will be necessary to obtain R 2 3 (among other items). 


Now 


^13 “ T 23 T n 


T \2 nxl 
T n mxl 
Tn kxl 
r 23 nxm 


F-2 


* Note these are “cumulative” voltage transfer rates from a common source point and are available (computed 
and saved during outbound computations). Using them instead of multiplicative chain multiplies could be a 
computational advantage during the reverse direction computation, albeit the “cumulative” ratios need to be 
computed during the outbound processing, but ti?ese can be done once and for all. 


MO 
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Obviously if m> n then T n ” [ does not exist. In this case, the above equation can be por- 
tioned as indicated by the following equation {as pointed out by J.T. Gajjar). 


13 hi hi 


nxm 

nx{l—m) 

BU 

nxm 

X 

mxm 

nixil— m) 

t hp 

T\3R 




t \ip 

Tur 


Then 


- T 22 

’A 2 P 

F-4 

r>3R - T 2} 

T 12 R 

F-5 


Now T [2p is square and presumably 7" 12 / r l exists so that 

hiR “ T\ 2 p T[ 2 p~ ] 

/ i 

Obviously by recording the rows and columns of T i 3 « T 22 T n could generate ^ t 

combinations for T\ 2P , he., the number of equations in T 13 — T 22 — T n represent an over- 
determined case for determining r 23 . In practice those columns (of T n ) associated with the 
m terminals of E 2 will be used. Note that ihe ability to make this identification is possible 
since a record is kept during the outbound processing of any transitions (3:2, 2:1, 3:1) and 
transpositions which may have been involved in the 7} 2 path. 
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USING RECIPROCITY TO COMPUTE REVERSE DIRECTION 
VOLTAGE TRANSFER RATIO MATRICES 

R.C. Rustay 


The following attachment is a verbatum copy of an internal memo concerning the relation 
ship of voltage transfer ratio matrices, driving point admittance matrices and the transadmit 
tance matrices. It is included in this report for possible tutorial interest. 


DIFNAP SYSTEM MEMO # 21 
10/14/80 


TO: C.A. Sum 
J. Gajjar 
W, Hughes 
R. Wooding 
R. Rustay 

FROM: R.C. Rustay 

SUBJECT: Comments on Voltage Transfer Ratios with Potential Application to the DIFNAP 
Main Program NTWKANS 

Recent discussion with W. Hughes regarding the mathematical properties and relationship 
between outbound and inbound voltate transfer ratios have stimulated the writing of this 
memo. Also, recent discussions with J. Gaijar regarding my desire to compute the outbound 
voltage transfer ratio matrix from the source to any section, also has stimulated comments in- 
cluded herein regarding the application of these analytical results to reverse path computa- 
tions. 
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« , , . l l M i 

Reciprocity of Transimpedance and Transadmittance 

The following is a short tutorial discussion of transimpedance and transadmittance recipro- 
city. This discussion will begin assumming reader acceptance of the completeness and gen- 
erality of the "y" parameter representation of linear circuits for those which "y" parameters ex- 
ist. Consider the following circuit: 



= YE 
mxm 


Figure 1 


where / is an m vector corresponding to the currents entering the m terminals. Similarly E is 
an m vector of the terminal voltages measured relative to a (zero potential) ground. Next 
separate the vectors / and £each into two pardoned subvectors / t , I 2 and E\, E 2 , i-e. 


I = 


h 

h 


E = 


Ex 

Ei 


where I x and E\ are n vectors associated with the “first” n terminals and simmilarly / 2 and 
E 2 are m-n vectors associated with the remaining m-n terminals. Then Figure 1 can be 
redrawn to 



n x n 


(m-n) x (m-n) 


Figure 2 
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where the Y mm matrix has been partioned as shown. T l2 r , the transpose of Y n occurs be- 
cause of the intrinsic symmetry of matrix K n , Y T =* Y, Figure 2 can be redrawn to frame it 
into a more conventional two port, i.e. an input and output 



Figure 3 


The vector matrix equation associated with Figures 2 and 3 are 



h 


I'll ^ 


l 2 =Y a %+Y 22 E 2 

h 


^12 T ^22 

X 


from which all subsequent results can be obtained. 


1 ) 

2 ) 


Trans admittance 

Suppose all the E 2 terminals are shorted to ground, i.e., £ 2 =0 then we observe from 
equation 1 that 7 n is the input driving point admittance* with the output terminals shorted 
and that the transadmittance* is Y\ 2 T (from equation 2). Similarity if input terminals are 
shorted to ground, i.e., 2^ = 0 then the driving point admittance* seen looking into the output 
terminals is T 22 and the (output to input) transadmittance* is seen to be F 12 . Hence the tran- 
sadmittances between the two parts are said to be reciprical, i.e. one is equal to the transpose 
of the other. 


Transimpedances 

The equations 1) and 2) can be inverted to the form 



*i 


%\] %12 


St-Zu T I\+Znh ** 

E 2 


^12 r %22 

X 


where Zn, Zi 2 , Z 22 are the corresponding matrix partitions of Zand 

Z - Y~ l 5) 


Since Y is symmetric, so also is Z. 

Now suppose all the E 2 terminals are open circuited, i.e., / 2 a “ 0 and a set of input currents 
are caused to occur. Then from equation 3) Z u is seen be the input driving point im- 

* all admittances are in general matrices which can be lx 1 
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pedance while from equation 4) 2’ 12 7 ' is seen to be the transimpedance (input to output). 
Similarly with the input terminal open circuited, i,c„ /j-0, then Z 22 is the driving point im- 
pedance looking into the output port and Z\ 2 is the (output to input) transimpedance, Hence 
as before the transimpedances are also seen to recipricol, i.e., one is equal to the transpose of 
the other. 

Note that the matrix elements in equations 3 and 4 are the equally well known “c" 
parameters which can obviously be expressed in terms of the "y" parameters, 

Relationship of ZP W Z n , Z n to Y n , Y n , Y n 

For possible future reference the relationship of Z n , Z l2 , Z 2V will now be given. These 
relationships can be obtained by matrix algebra operations on equations 1) and 2) or 
equivalently making use of well known partioned matrix inversion procedure, the following 
results can be obtained 

Zn- «) 

Z n --[ I'll- y a Zn' I'i2 7 ]" 1 r\i Zn [ - - Zn ^2 1 * 

[I'H-l'i^l'if'l'ia]- 1 *> 


Admittance Matrix Equivalent for a Two Port Circuit With External Loading 

The purpose of discussion is to show that a two port linear circuit, with external loading, 
i.e., 





E 


2 


Figure 4 

can be easily converted to an equivalent unloaded two port of the type discussed in the 
preceeding section. Writing obvious equations 

/, - r, £,+/,* - (r,j+ /,)£,+ y 12 e 2 9) 

h- Y 2 E 2 +I 2 >- Y u r E l +(Y 22 +Y 2 )E 2 10) 

Hence, the equivalent K matrix is seen to be 


* Many other equivalent expressions can be formulated 
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V s i ; 

' *12 

I 



Y ^ 
X 12 

Y + Y 
| x 22 r *2 






A 

E 


Figure 5 


Thus the discussions of the proceeding and following sections are perfectly general. 

Relationship Between Input to Output and Output to Input Voltage Transfer Matrices 

This section is concerned with the properties of the voltage transfer ratio matrices from in- 
put to output and output to input for a two port linear circuit of the configuration shown in 
Figure 3 where it is assumed that any external loading has been imbedded into the network Y 
matrix as discussed in the preceeding section. This discussion also assumes observability and 
controlability (also implied for the preceeding sections but more likely to be overlooked in 
this context) which imply that r n _1 and Yif X exist. 

Voltage Transfer Matrix G n from Input to Output Port 

Given a voltage E x (vector) applied somehow to the input port of Figure 3, the voltage E 2 
at the output port with those terminals open circuit, i.e. / 2 «-0 (remember any loading on the 
output port has been imbedded into the circuit Y matrix*) is easily seen to be from equa- 
tion 2) 

£ 2 - Y n -'Y n TE x ~ G n E x I 2 - 0 11) 

i.e., 

g, 2 - - r 22 -> r ia r 12 ) 

Alternatively using the "Z” representation of equations 3) and 4) it is easy to obtain, for the 
same condition that / 2 -0 

E 2 -Z n T Z n - x / 2 -0 13) 

i.e., 

G n -Z n T Z x c' 14) 


* i.e., from the preceeding section this output loading admittance would be linearly added to the original circuit 
and thus Y u ~ x in equations 11) and 12) will involve this loading 
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Voltage Transfer Matrix <7 2 i from Output to Input Port 

Similarly, given a voltage E 2 (vector) applied somehow to the output port of Figure 3, the 
voltage E\ at the input port with those terminals open circuit, i,e„ / j«0 (remember any load- 
ing on the input port has to be considered and also imbedded into the circuit Y matrix) is 
easily seen to be from equation 1) 

E x - Y n -' Y n E 2 - G n E 2 7 r 0 15) 

i.e„ 

G 2 1 - - Y u - ■ T,J 16) 

Similarly and alternatively using the “Z” representation of equations 3) and 4) it is easy to 
obtain for the same condition that /j-0 

E\ — Z\ 2 Z 22 ~ ^E 2 /j—0 17) 

i.e,, 

G 2 \ “ Z\ 2 Z 2 2~'^ IS) 

Observations Regarding Voltage Transfer Ratio Matrices 

Examining the voltage transfer ratio matrices, i.e., 

C|2-- Y 22 - l Y n T= Z l2 TZ { r 1 19) 

On-- Y n -'Y n - Z l2 Z 22 -' 20) 

observe that they involve the “transadmittance” F 12 or the “transimpedance” Z 12 commonly 
between them- Obviously G n and G 2l are not reciprocol. Also if G l2 , Y 22 , and Y n (or G n , 
Z 22 , Zii) are knov/n then (7 21 can be calculated, i.e., 

C 2I - Yn'G n TY 22 = Z n G n TZ n -' 21) 

Similarly for G\ 2 (but not usually of interest) 

<3,2 - Y 22 -'G 2 JY n => Z 22 G 2 JZ n -' 22) 

It is appropriate at this point to reemphasize that as used in these discussions concerning 
voltage transfer ratio matrices, F u , Y 22 and the Z (J contain the embedded effects of any ter- 
mination (see Figures 4 and 5) (including any “source” admittance) since obviously the <7 12 
and G 2 i are dependent, in the context implied here, on both terminations. 

The alternate formulations given above, in terms of For Z parameters represent short cir- 
cuit (admittance) and open circuit (impedance) approaches respectively. The latter form is 
used in the DIFNAP programs because open circuit conditions (with terminations imbedded) 
are generally of interest, i.e. terminal voltages. In the context and terminology of the DIF- 
NAP system Z u -1 is referred to as Y,„ and Z 22 ~ l is referred to as Y out so that equations 21) 
and 22) become 

G 2 i- Yin~ l G\ 2 T Y 0Ut 23) 


G 


12 


Yout~ x G 2 jY ln 


24) 
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It is interesting to observe the physical interpretation of Z form of Equations 19) and 20). 
Consider the voltage transfer C? l2 from port 1 to port 2. The current /j into port l is 
/i - Y i „E l and the port 2 voltage E 2 due to /, is £ 2 - Z 12 % so that E 2 - Z n T Y m E\ or 
G, 2 - Z l2 TY m as before. 


The relationships 21 and 22 involving Z\\ - Y tn ~ [ and Z 22 - Y ou ~ x offer a procedure for 
calculating the reverse path voltage transfer ratio matrices in the DIFNAP program 
NTWKANSi. Some discussions on such an algorithmic procedure are discussed in the next 
section. 


Comments Regarding Algorithmic Procedure for Computing Reverse Path Voltage Transfer 
Matrix 


This discussion is concerned with some tentative and embryonic ideas regarding algorithm- 
ic procedures which could be implemented into the DIFNAP computer program NTWKANSi 
for computing the reverse path voltage transfer ratio. 

Consider, for the purpose of an example to illustrate the possible algorithm, the following 
cascaded sections. (It is obvious that branching presents no conceptual difficulties; practically 
it will depend upon having available appropriate matrix admittances associated with 
branching.) 



As a result of normal outbound response calculations, the voltage transfer ratios G 12 , G 23 , 
G34, the network driving point admittance K //V and the “equivalent Norton source admit- 
tances” Y out , Yq UT i: , Y 0( j Ti are available. (As mentioned above there may be a lack of com- 
plete saved data when branching is involved.) Now consider the three following cases: 

<?21 

From equation 23 

^ 21 “ G\ 2 t Y 0 ut { 25 ) 

^31 

Similarly 

G 3\ = Yin~ V Gu t Yovt 2 26) 

But 

G u - G 23 G u ~> G l3 r =* G n TQ 23 71 27) 


^31“ ^~ l( ^12 r< ^23 r ^Ot/r 2 

G 4\ 

Similarly 


28) 
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^ 41 * Yin“ 1g \4 T YoUT 3 
But 


29 ) 


<?14 - C7 34 C7 23 C?i 2 -> C?,4 r - 30) 


^41- YiH^Gn^G^Your, 31) 

Hence, it is easily observed that a recursive computational algorithm can be realized in the 
same sense that the present Norton’s equivalent source admittances (looking in the reverse 
direction) such as Y 0UTx , Y 0UTl , and Y 0UT} represent the total affect of the network upstream 
from that point towards the source. For example during the normal outbound processing, say 
to get the voltage *? 3 , the result of the computation 

hvr'Gn T Gv T 32) 

could be saved with the other data saved with section 2. Then progressing to section 3 the 
corresponding quantity 

y W r' G n T (!2 } T 0 3t r~ <.Y INi -'G n TG 21 T) * G } / 33) 

could be computed and saved with that section’s data. (Note an appropriate RCR subroutine 
already exists for computing C-AB T .) 

Then either on demand, or at the time the above computations are made (and saved for 
later recall) the reverse path voltage transfer ratio could be computed, i.e., for example 

G41 " (J7/v" : I< 3 ! 12 T G23 T G 32 0 * Y0UT1 34) 

v "tusly many “bookkeeping” activities need to be considered. For example, during the 
no,„ abound NTWKANSi calculation, the recursively accumulated matrix such as 
Y/ff l Gi 2 T G 2 2 T should be saved in core and saved with the data associated with the next 
downstream section (or two with branching) to avoid an extra disc access when later recalling 
this data for computing on demand the reverse direction voltage transfer ratio matrix. How- 
ever, if the extra computing time is acceptable, the reverse direction voltage transfer ratio 
matrix could be computed for every section. However, it is not necessary that all sections 
will be evaluated tending to “throw out” this procedure. For the “on demand” procedure an 
extra dimension would be required to the LIFO (push down stack) matrix array to accommo- 
date to the potential level of branching traversed while computing an outbound propagation. 

Finally it is noted that the above procedure implies transfers back to the original source. 
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APPENDIX 6 

MAIN PROGRAM NETGENS1 FOR GENERATING GENERIC NETWORKS 

R,C. Rustay 

Included in this Appendix is a copy of an internal memo describing a main program NET- 
GENS 1 for generating generic networks 


DIFNAP MEMO #24 

n /i 7/80 


Memo To 


From: 

Subject: 


C, 

J. 

R. 

0 . 

W„ 

R. 

R. 


A. Stutt 

Gajjar 

Wooding 

Rodriguez 

Hughes 

Rankin 

Rustay 


R. C. Rustay 

Main Program NETGENSI for Generating Generic Networks 


Please be advised that a main program NETGENSI is now opera- 
tional. This program provides a convenient procedure for generating 
a simple class of generic network files, of the form DNWKINij, and 
which are characterized by the following items: 


1) The generated network is everywhere three phase 

2) Only one level of branching is available, ie. as illus- 
trated on the attached sketch. 

3) The networks shown in the sketch are composed of segments, 
the one originating at the source or root section being 
designated as the main "route" segment; all others being 
designated as laterals with up to a maximum of 10 being 
user definable. 

4) Each segment will have an invariant line/cable configura- 
tion which is user definable. 

5) Each segment can be composed of many sections with the 
sections in any segment being equal length (except possibly 
for the last one). These section lengths are user definable 
for the main route and all laterals. 

6) The user can invoke an option to cause the program to auto- 
matically place a distribution transformer at the end of 
each section. The program automatically balances this 
loading by cycling the assignment of the DT in turn to each 
phase. The user can elect to defeat this option in which 
cose all lines will be unloaded. 
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The user defines the total length of all segments. 

The user defines the distance along the main route (from 
the root/source) at which each lateral will be attached. 
Note ; In the present first edition, the logic necessary 
to introduce "zero length" sections to permit tertiary or 
higher order branching has not been included. Therefore, 
at the present time, it will be the user's responsibil ity 
to preplan to avoid this situation. If a need occurs and 
time/resources allow, this additional logic for higher 
order branching could be Included. 

The user can arrange for a lateral to be appended to the 
end of the original main route segment by specifying 
connection distance for a lateral to be greater than the 
total length of the main route. This procedure can be 
cascaded. This procedure allows the user to construct 
a cascaded sequence of segments where each segment may 
have Its own and different line types and/or section 
lengths. Subsequent laterals can be specified to branch 
off anywhere on the, now, augmented main route. 

The resulting network is placed in a file whose mane is 
of the form DNWKINij where ij is specified by the user. 
The program includes a check/alert to the user to guard 
against over-writing an existing file of the same name. 

The output file includes an additional column which gives 
the accumulated length along any segment (including any 
augmented main route). This added column will not inter- 
fere with the reading of this file by the main network 
program NTWISERSj; this last column is simply ignored. 
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APPENDIX 7 

MAIN PROGRAM SUBNETS1 
R.C. Rustay 

Included in this Appendix is a copy of an internal memo describing a main program SUB- 
NETS 1 which can extract as a subnetwork a portion of a larger network. 
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DIFNAP SYSTEM MEMO #23 


11/13/80 


Memo To: C. A. Stutt 

J. Gajjar 
W. Hughes 
R. Rankin 
D. Rodriguez 
R. Rustay 
R. Wooding 

From: R. C. Rustay 

Subject: Availability of Main Program SUBNETSI to Extract Subnetworks 

Please be advised that a main program SUBNETSI is now opera- 
tional. This program allows the user to extract from a given network 
file (usually with a file name like DNWKINij) a "subtree" of that 
network. This subtree is defined to the program by specifying the 
section number (which, of course, must already exist) which is to 
become the new root section. The program then extracts the portion 
of the original network which is outbound from that new root section, 
and places the resulting subnetwork specification into a file whose 
name is specified by the user. The program ‘SUBNETSI employs conversa- 
tional input and is self-explaining during running of the program. 

Conforming to RCR policy, the user can and should exit from the 
program by typing a negative number when prompted. 

Any section identification number cannot exceed 2000 with the 
present dimensioning. This arbitrary limit can easily be increased, 
if necessary, at the expense of additional core requirement. 


RCR: jw 
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APPENDIX 8 

NETWORK DESCRIPTION 

This Appendix contains specific information on the parameter variations analyzed during 
the study. Each variation is identified by its network (or system) number which is also used 
to identify the output data for that variation. Unless a parameter is specifically mentioned, it 
can be assumed that the parameter is the same as it was for the associated reference network. 
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Reference Networks 


Network #8 


Network #9 


- 1/4 Wavelength underground cable network. 

. Same as Network #0 except section lengths 
reduced to 42.5m for a total length of 
3400m. 

- Nominal underground cable network 

. Aluminum phase conductors GMR - .385 inches 

. 21 copper neutral strands, each with a radius 

of .64 inches, evenly spaced around the phase 
conductor at a radius of .896 inches. 

. Main insulation has an inner radius of .515 
inches and an outer radius of .750 inches. 

. Outer jacket has an inner radius of .780 inches 
and an outer radius of .880 inches. 

. Center to center spacing of phase cables is 
2.0 inches, equilateral spacing assumed. 

. Phase conductor temperature - 75 e C 

. Neutral temperature = 45°C 

. Earth resistivity = 100 n / m 3 

. Signal frequency =8.13 kHz 

. Total of eighty transformers, one loading each 
section, rated at 25 kva with load admittances 
of -78.8 dB at -27.3 degrees. 

. Eighty sections, each 200m long for a total 
length of 16,000m. No branches, no loops, 
entire line is three phase. 

. Diagonal load at the end of the line has three 

-4 

equal elements, 4.77 x 10 mhos -68.22° each. 
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Network #10 


Nominal overhead reference network. 


. Geometric arrangement as shown in Fig. (6a). 

Phase conductors are 336 Widgeon ACSR with 1/0 

ALCOA ACSR neutral wire. Diameters are .721 inches 

and .398 inches, respectively. Cross-sectional 

2 2 

areas are .408 inches and .124 inches , respectively. 
This line type code is 101. 

. Total of eighty transformers, one loading each 
section, rated at 25 kva, with load admittances 
of -78.8 dB at -27.3 degrees. 

. Eighty sections, each 200m long, for a total 
length of 16,000m. No branches, no loops, entire 
line is three phase. 

. Temperature of conductors = 25°C 

. Earth resistivity = 100 fi / m 3 

. Signal frequency =8.13 kHz 

. Diagonal load at the end of the line has three 
equal elements, 4.77 x 10“^ -68.22° mhos each. 


8.2 Overhead Geometry Variation Networks 


System #11 


System #12 


System #13 


Referred to network #10; Phase conductor changed 
to 4/0 ALCOA ACSR. Diameter reduced from .721 
inches to .563 inches, a 21.9% reduction. Cross 
sectional area reduced 39.0%. 

Referred to network #10; Phase conductor changed 
to 2/0 ALCOA ACSR. Diameter reduced from .721 
inches to .447 inches - a 38% reduction. Cross- 
sectional area reduced 61.6%. 

Referred to network #10; Phase conductor changed 
to 447 HEN ACSR. Diameter increased from .721 
inches to .883 inches, a 22.5% increase. Cross- 
sectional area increased 50.0%. 


8-3 




9W l 


System #14 

- Referred to Network #10; Neutral conductor changed 
to #4 ALCOA ACSR. Diameter reduced from .398 
inches to ,250 inches, a 37.2% difference. Cross- 
sectional area reduced 60.5%, 

System #15 

- Referred to network #10; Phase conductor to phase 
conductor spacing increased from 3.667 feet to 
4.033 feet, a 10.0% difference. 

System #16 

- Referred to network #10; Phase conductor to earth 
spacing increased from 35.063 feet (for Phase A) 
and 33.667 feet (for phase B & C) to 38.569 feet 
and 37.033 feet, respectively. This is an increase 
of 10.0%. 

System #17 

- Referred to network #10; Neutral* conductor to earth 
spacing increased from 27.333 feet to 30.066 feet, 
a 10.0% difference. 

System #18 

- Referred to network #10; Miscellaneous conductor 
added, 23.458 feet above the earth. Assumed to be 
at ground potential . 

System #19 

- Referred to network #10; Two miscellaneous conductors 
added at 23.208 and 23.458 feet above the earth. 

Both assumed to be at ground potential. 

System #20 

- Referred to network #10; Two open circuit secondary 
conductors added at 26.000 and 26.667 feet above the 
earth. Both assumed to be at ground potential. 

System #21 

- Referred to network #10; Neutral conductor changed 
to Triplex. Diameter increased from .398 inches to 
.92 inches, a 131% difference. Cross-sectional area 
increased 434%. 
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System #22 
Network #28 
Network #29 
Network #30 
Network #31 
Network #32 - 
Network #33 - 
Network #34 - 
Network #35 - 
Network #36 - 
Network #37 - 
Network #38 - 
Network #39 - 


- Referred to network #10; Configuration changed to 
spacer cable. See Figure (M) 

- Referred to 1/4 wavelength version of network #10; 
same variation as system #11. 

- Referred to 1/4 wavelength version of network #10; 
same variation as system #12. 

- Referred to 1/4 wavelength version of network #10; 
same variation as system #13. 

- Referred to 1/4 wavelength version of network #10; 
Same variation as system #T4. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #15. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #’6. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #17. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #18. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #19. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #20. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #21. 

Referred to 1/4 wavelength version of network #10; 
same variation as system #22. 
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Figure (8-1) 
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8,3 Transformer Loading Variation Networks 


Network #13 - Referred to network #10; The diagonal load at the end 
of the line is removed, leaving it open circuited. 

Network #14 - Referred to network #10; The admittance of each 

distribution transformer load is doubled to >72.8 dB 
at -27.3 degrees . 

Network #17 - Referred to network #10; The diagonal load at the end 
of the line is removed while at the same time the 
admittance of each distribution transformer load is 
doubled to -72.8 dB at -27.3 degrees. 

Network #42 - Referred to 1/4 wavelength version of network #10; 
same variation as network #13. 


Network #45 - Referred to 1/4 wavelength version of network #10; 
same variation as network #17. 

Network #56 - Referred to 1/4 wavelength version of network #10; 
same variation as network #14, 


Network #57 


Referred to 1/4 wavelength version of network #10; 
all distribution transformer loading is neglected. 


Network #58 


Referred to network #10; same variation as network #57. 


Network #59 - Referred to network #10; admittance of each distribu- 
tion transformer load is cut in half to -84,8 dB at 
-27.3 degrees . 

Network #60 - Referred to network #10; overhead geometry varied as 
in system #11, loading varied as in network #59. 

Network #61 - Referred to network #10; overhead geometry varied as 
in system #12, loading varied as in network #59. 
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Network 162 - Referred to network #10; overhead geometry varied as in 
system #13, loading varied as in network #59. 

Network #63 - Referred to network #10; overhead geometry varied as 
in system #14, loading varied as in network #59. 

Network #64 - Referred to network #10; overhead geometry varied as 
in system #15, loading varied as in network #59. 

Network #65 - Referred to network #10; overhead geometry varied as 
in system #16, loading varied as in network #59. 

Network #66 - Referred to network #1 Op overhead geometry varied as 
in system #17, loading varied as in network #59. 

Network #67 - Referred to network #10; overhead geometry varied as 
in system #18, loading varied as in network #59. 

Network #68 - Referred to network #10; overhead geometry varied as 
in system #19, loading varied as in network #59. 

Network *9 - Referred to network #10; overhead geometry varied as in 
system #20, loading varied as in network #59. 

Network #70 - Referred to network #10; overhead geometry varied as 
in system #21, loading varied as in network #59. 

Network #71 - Referred to network #10; overhead geometry varied as 
in system #22, loading varied as in network #59. 

8.4 "Set up" Time Reduction Networks 

Network #11 - Referred to network #10; the length of every other 

section is doubled to 400m with the in-between sections 

reduced to zero length. Same overall length and loading 

are maintained. Number of sections is effectually cut 

in half, though not actually. [ ' 

r •* 

1 ; *• 

t 


r 

f. 

ta 
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Network # 12 - Referred to network #10; the length of every third 

section is tripled to 600m with the in-between sections 
reduced to zero length. Because 80 is not evenly divi- 
sable by 3, the next to last section is doubled instead 
of tripled in length to maintain the same overall length 
as the reference network. Same loading is maintained. 
Number of sections is effectively reduced to almost a 
third of the original number. 

Network #15 - Referred to network 10; same variation as network #11 
with the admittance of each transformer load doubled 
to -72.8 dB at -27.3 degrees. 

Network #16 - Referred to network #10; same variation as network #12 
with the admittance of each transformer load doubled 
to -72.8 dB at -27,3 degrees. 

Network #18 - Referred to network #10; Temperature changed from 35°C 
to 30°C, a 20.0% increase. 

Network #19 - Referred to network #10; Temperature changed from 25°C 
to 20°C, a 20.0% decrease. 

Network #20 - Referred to network #10; temperature changed from 25°C 
to 10°C, a 60.0% decrease. 

Network #21 - Referred to network #10; Temperature changed from 25°C 
to 0°C, a 100.0% decrease 

Network #22 - Referred to network #10; temperature changed from 25°C 
to -10°C, a 140.0% decrease. 

Network #23 • Referred to network #10. Earth resistivity changed 
from 100 n / m 3 to 50 n /, n 3, a 50.0% reduction. 

Network #24 - Referred to network #10. Earth resistivity changed 
from 100 n / m 3 to 75 n / m 3, a 25.0% reduction. 
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Network #25 - Referred to network #10. Earth resistivity changed from 
100 '7 m 3 to 125 ‘7 m 3, a 25.0* Inemse, 

Network *26 - Referred to network #10, Earth resistivity changed 
from 100 l! /|||3 to 15Q '7 ra 3. a 50.0* Increase. 

Network #27 - Referred to network #10. It is not assumed that the 
neutral conductor is at ground potential everywhere 
along the line, 

Network #40 * Referred to 1/4 wavelength version of network #10, 

Same variation as network #11 except that section 
lengths are adjusted for 1/4 wavelength. 

Network #41 - Referred to 1/4 wavelength version of network #10. 

Same variation as network #12 except that section 
lengths are adjusted for 1/4 wavelength, 

Network #43 - Referred to 1/4 wavelength version of network #10. 

Same variation as network #15 except that section 
lengths are adjusted for 1/4 wavelength, 

Network #44 ~ Referred to 1/4 wavelength version of network #10, 

Same variation as network #16 except that section 
lengths are adjusted for 1/4 wavelength. 

Network #46 - Referred to 1/4 wavelength version of network #10. 

Same variation as network #18, 

Network #47 - Referred to 1/4 wavelength version of network #10, 

Same variation as network #19, 

Network #48 - Referred to 1/4 wavelength version of network #10, 

Same variation as network #20. 

Network #49 * Referred to 1/4 wavelength version of network #10, 

Same variation as network #21, 


t 

ll 

I 

ii 

1 

Network #50 - Referred to 1/4 wavelength version of network #10. } 

Same variation as network #22. I 

Network #51 - Referred to 1/4 wavelength version of network #10. 

Same variation as network #23. 

Network #52 - Referred to 1/4 wavelength version of network #10. 

Same variation as network #24. 

Network #53 - Referred to 1/4 wavelength version of network #10. 

Seme variation as network #25. 

Network #54 - Referred to 1/4 wavelength version of network #10, 

Same variation as network #26. 

Network #55 - Referred to 1/4 wavelength version of network #10. 

Same variation as network #27. 

Network #74 - Referred to network #10. Each section is reduced 
in length to 190m to give an overall length reduc- 
tion of 800m, or 5%. 

Network #75 - Referred to network #10. Each section is increased 
in length to 210m to give an overall length increase 
of 800m, or 5%, 

8.5 Cable Networks 

Network #76 - Referred to network #9. Earth resistivity changed 
from 100 ^/ m 3 to 50 n / m 3, a 50.0% reduction. 

Network #77 - Referred to network #9. Earth resistivity changed 
from 100 n / m 3 to 150 n / m 3 , a 50.0% increase. 

Network #78 - Referred to network #9. Earth resistivity changed 
from 10Q n / m 3 to 200 n / m 3, a 100.0% increase. 

Network #79 - Referred to network #9. Temperature of phase 

conductors changed from 75°C to 65°C, a 13.3% decrease. 

Temperature of neutral strands changed from 45°C to 

35 °C, a 22.2% reduction. 
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Network #81 - Referred to network #9. Temperature of phase 

conductors changed from 75°C to 95°C, a 26.6" increase. 
Temperature of neutral strands changed from 45°C to 
65°C, a 44.4% increase. 

Network #82 - Referred to network #9. The length of every other 
section is doubled to 400m with the in-between 
sections reduced to zero length. Same overall length 
and loading are maintained. Number of sections is 
effectively cut in half. 

Network #83 - Referred to network #9. The length of every third 
section is tripled to 600m with the in-between 
sections reduced to zero length. Because 80 is not 
evenly divisible by 3, the next to last section is 
doubled instead of tripled in length to maintain the 
same overall length as the reference network. Same 
loading is maintained. Number of sections is effec- 
tively reduced to almost a third of the original number. 

Network #84 - Referred to network #9. The admittance of each dis- 
tribution transformer load is doubled to -72.8 dB 
at -27.3 degrees . 

Network #85 - Referred to network #9. The admittance of each dis- 
tribution transformer load is cut in half to -84.8 dB 
at -27.3 degrees . 

Network #86 - Referred to network #9. All distribution transformer 
loading is neglected. 

Network #87 - Referred to network #8. Same variation as network #76. 

Network #88 - Referred to network #8. Same variation as network #77. 

Network #89 - Referred to network #8. Same variation as network #78. 
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Network 

#90 - 

Referred 

to 

network 

#8. 

Same 

variation 

as 

network 

#79. 

Network 

#91 - 

Referred 

to 

network 

#8. 

Same 

vari ation 

as 

network 

#80. 

Network 

#92 - 

Referred 

to 

network 

#8. 

Same 

variation 

as 

network 

#81 . 

Network 

#93 - 

Referred 

to 

network 

#8. 

Same 

variation 

as 

network 

#82 



except that 

section 

lengths are adjusted 

for 1/4 wave- 


length network. 

Network #94 - Referred to network #8. Same variation as network #83 
except that section lengths are adjusted for 1/4 wave- 
length network. 


Network #95 
Network #96 
Network #97 
Network #98 


Referred to network #8. 

Referred to network #8, 

Referred to network #8. 

Referred to network #9. 
reduced to 190m to give 
of 800m, or 5.0%. 


Same variation as network #84. 

Same variation as network -#85. 

Same variation as network #86. 

The length of each section is 
an overall length reduction 


Network #99 - Referred to network #9. The length of each section is 
increased to 210m to give an overall length increase 
of 800m, or 5 .0% . 
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APPENDIX 3 

Y PARAMETER ANALYSIS OF SYMMETRIC DISTRIBUTION TRANSFORMER 
WITH BALANCED LOADING AND BRIEF DISCUSSION 
OF A RLC LUMPED PARAMETER MODEL 

R.C. Rustay 


The material contained in this Appendix was originally issued as an internal memo early in 
the Phase II work as an attempt to summarize, at that time, what in our judgement was apt to 
be the most useful type of distribution transformer model (and also any single phase two 
winding transformer). At that time, two candidate models; 

Y Parameter Derived Algebraic Model with 
Frequency Dependent Parameters 

A RLC Lumped Parameter Model 

were selected. 

The writing of the original memo followed laboratory measurements and analytical studies 
performed by J.T. Gajjar which demonstrated the feasibility of the Y parameter derived alge- 
braic model. Also, at that time, R.W. Rankin had performed response measurements (see 
plot at the end of this Appendix) which strongly suggested the feasibility of a simplified RLC 
lumped parameter model for frequencies less than, say, 50-100 kHz. R. Wooding subse- 
quently showed that indeed an RLC model, with appropriate selection of the lumped parame- 
ter values could account for the measured response. Subsequent “full scale” measurements 
were made to develop the parameters for both models. Also at that time, measurements 
were made of critical admittance and transfer characteristics which were used to verify the 
selected model. 
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DIPNAP SYSTEM MEMO #14 


MEMO TO: C.A. Stull 

J. Oajjur 
R. Woudiug 
R, Rankin 
W, Hughes 

FROM: R.C. Rustay 

« 

SUBJECT: 3' Parameter Analysis of Symmetric Distribution Transformer 

with Balanced Loading and Brief Discussion of a RLC 
Lumped Parameter Model 

1. INTRODUCTION 

This memo presents, for reference, the results of U Y” parameter analysis of a distribution 

transformer for the special conditions of symmetric construction and balanced secondary loading. 

These conditions were suggested by J. Gajjar for modeling/simulation because: 

1. in most cases the algebra is greatly simplified and leads to inverse linear relationships involv- 
ing loading 

2. these conditions/assumptions are the best that can be assumed for modeling in the sense that 
there is no apriori basis for assuming other than balanced secondary loading and that for the 
lower frequencies at least the DT’s do appear to be electrically symmetric with respect to the 
secondary electrical characteristics 

3. these conditions and the resulting simple algebraic models reveal and suggest easy to measure 
sets of parameters (open circuit voltage transfer ratios, open circuit impedances) us well as 
some of the conventional Y parameters which avoid the sensitivity to measurement 
precision/error that appear to occur when all "v" parameters are used. 

4. these conditions lead to particularly convenient and computationally stable algebraic formula- 
tions for driving point admittances 

Also presented are comments on a “zeroeth” order two winding RLC (plus ideal transformer model). 
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2. ALGEBRAIC MODEL DERIVED FROM “ F” PARAMETER ANALYSIS 

The muierial contained in this memo represents a simplification of the more general results con- 
tained in Appendix Q of Vol. 3 of the “RF Model — ' " final report, plus rearrangements to include 
other than V Y" parameters. 

Figure l-l is the generic/gencral distribution transformer schematic model with connected loads 
from which the various special subcases can be derived. This schematic serves to define terminology 
and convection labeling. 

Y | Admittance load connected across primary terminals 

Y}, Y j, Y 2 j Admittance representation of secondary load 

y iV Neutral to ground admittance 



Figure 1 - 1 . General Distribution Transformer Schematic Model with Connected Loads 


Figure 1-2 is a general distribution transformer schematic (without loading) and serves to define its 
associated y tJ parameter representation in terms of the associated Y t admittance matrix shown on 
Figure 1-3. In Figure 1-3 only the independent matrix elements y tj are shown; the other matrix ele- 
ments can be derived/deduced from symmetry and the (row-column) zero sum property of an 
ungrounded network. Figure 1-4 is similar schematic for the terminal loads, including the earth to 
ground admittance Y v . Figure 1-5 is the associated terminal load matrix Y L . The matrix admittance 
representation corresponding to Figure 1-1 is then seen to be the sum Y L +Y/-. 

Various special cases can be easily obtained using the admittance matrix Y L + Y/- and applying the 
terminal conditions associated with the special cases. 

Finally, it should be noted that the convention adhered to in this memo has all terminal voltages e, 
representing terminal to ground potential, and all terminal currents /, positive into circuit as shosvn on 
various schematics. ; 

i 

As mentioned above, this memo is concerned with results associated with electrically symmetric r : 

transformer characteristics and balanced secondary loading. These conditions are represented by the f 

statements: jl 


3-3 





ORIGINAL PA® K 
OF POOR QUALITY 



Figure 1-2. General Distribution Transformer Schematic 


y 12 

y t3 

y 14 

3*22 

^23 

^24 

-- 

^33 

>>34 

-- 

— 

^44 


Figure 1-3. Associated Distribution Transformer Admittance Matrix Yp 

7j- Y 2 
Yu - Y 22 
Yu = -Ti2 
^44 = Y u 
Y 2 4 - Yu - -Y\ 2 
y 34 = r l2 

leaving only the following “independent” elements 

yn yi 2 yn yu ^14 

Figure 1-6 shows the associated distribution transformer admittance matrix Yp subject to these electrical 
symmetry conditions. 
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Figure 1-4. Terminal Load Schematic 
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Y 2 +Y 23 
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-y 2 
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0 
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-Yl 

0 

0 

Yi 
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0 

-y 2 
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0 

Y |-f Y 2 + j 


Figure 1-5. Associated Terminal Load Admittance Matrix Y L 
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Figure 1-6. Associated Distribution Transformer Admittance Matrix Yr Symmetric Llectrical Properties 
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Case A - Phase to Neutral Primary Connection* N mural at Zero Potential 

This ease represents the most frequently occur my connection and assumption regarding the neutral 
being at zero potential. This case can be derived Irom die general case ot Figure 1*1 by applying the 
terminal conditions * 


£4 *“ **3 m 0 t Y\ — oo ) 


The resulting schematic is shown on Figure A*1 and the associated total admittance matrix V is gi\en 
in the following equation 


or in equation form 


>'tl+ Ft yu ~y\i 


Y « 


y ii y>>'b ) ;+ F;j ,Vxv~ Yu 

— 1 2 ,t*a*“^25 y’2i*F Yyr )‘>j 


ft 


JTrPl'i 3' tj 2 


h 

m 

yqj y 2:4* )YF jpjy— Fjj 

X 

h 


~*J1’|2 y 23"“ Yii Pj4* Fjj 



A* l 


A*2,A*3,A»4 



Figure A-l. Distribution Transformer Schematic - Phase to Neutral Primary Connection, Neutral at 
Zero Potential ( ir\— «») 

In accordance with the assumption stated earlier regarding a balanced secondary loading, it will be 
assumed for the remainder of this note that 


*3" 


With this balanced load assumption the equations become 


A-5 


• Which elVccnvcIy deletes die 4'th and S'tli row «u\d culutnus. 
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i\ 


y\&y i 

y i 2 

*"^12 


h 

- 

>'l2 

>'22+ ^2+ ^23 

>'23~ ^ 23 

X 

h 


~y 12 

-V23~ ^23 

>'22+ ^2+ Yl.\ 



By applying appropriate terminal conditions to these equations, all desired results can be obtained. 
Finally, it is evident that with symmetric electrical characteristics and balanced secondary loading, 
ej A ~e 2 Tor primary to secondary considerations.* Also it is of interest to note that for /| A 0 

f i - - (ej-ej) /i 4 0 A- 9 

>Ml+ M 

which is valid for any condition where both e 2 and e 2 (or their difference e 2 —e}) is specified by external 
conditions. 


* i.c., ij—i j=*0 
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Case A - Primary to Secondary: /’2“O“0> Yy* Y 2 , V |“0, Y\ -7 00 

This case represents the situation looking into the primary terminal 1 with Ki-0 (since we are 
interested in the input driving point admittance to the transformer). Applying the terminal conditions 
/ 2 WJ - 0 and the fact that ey - <*2 obtain easily 


— ey**e 2 m 


i l 


-^12 

>'22~>'23'h >^2+2 Yu 




.i'll 


2^2 


^22“^23+ ^2+2 Y 2 ) 


<?l 


If the open (secondary) circuit primary driving point impedance Z\\ 

J_ - * 

*11 


^11 


y,-o- y 


23 


y 22 y 23 

is available, J.T. Gajjar has suggested the following equivalent expression 


_L + 

L __l] 

^2+2/23 

z li 

y u ~ — 

Z U 

J'22""J'2?”t* )'2+2 Y 2 j 


e\ 


A- 10 


A* 11 


A- 12 


A-13 


Noting that the open (secondary) primary to secondary line to line voltage transfer ratio £12 is given by 

e 2 -e 2 -2y l2 


§ 12 




y22~~y2i 


Y 


23 


A- 14 


so that 


e 2 » — e 3 


2 Y 2 +2Y 2i 

S 12 ^12 


A-15 


' 1 




>'2+2 K 23 


^ 2+2 r 2 3- 


2^12 
£ 12 




A-l(> 


If the parameters z\\,yn,y\ 2 .g [2 are measured, the following identity-redundancy can be used to validate 
the measurements. 

— - yn + >'125-12 A- 17 

211 


It is seen that the dependancy on the secondary load condition Y 2 +2 Y 2 3 (remember Y$=* Y 2 ) has a 
simple inverse linear fashion. The forms A-15 and A-16 involve easily measured quantities 

g\ 2 ^t 2 * The form A-16 is intuitively and computationally attractive because — and y\[ 

’ ’ y r ' 1 

represent the asymtotic limits. Also further analysis will show that Yi+2Y 22 always has a positive 

£12 


real part and therefore will not ever be zero. 


* An equivalent and alternate expression is 


Y \\ + 

\- L--v 


1 11 

~ — yw 

-U 

1 - — 
y 12 

Y >+2 Yn 



2 


A-! 7a 
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Casft A - Secondary to Primary: /i*«0, Y^Y^ Y^—* 00 

The purpose of this condition is to reemphasize a conclusion reached earlier [see (A-9) that for any 
condition where e± and t*j or their difference e 2 -ej] is specified, then 

*1 -- 7^57 <*«-**> A. 18 

Because of various manners in which the secondary can be driven various special subcases will be con- 
sidered. 
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Case A - Secondary to Primary 1 1 — 0, f'yNjr* 0, Kj- )%—♦<» 

This case corresponds to u “floating"* drive with a voltage t’j-ei or where a ‘'balanced" drive with 
i?) A ~cj. In the former case obviously /j-wr'ir In the latter case 0 because of the tssumptions 
that Y- r* }S and the transformer us electrically symmetric. Then (as before) 




^12 

y\&Y\ 




A* 19 


h ^2 ^ . ^2 . T22-\‘\\l J’ 1*2 

‘ — M — ■ MM ) 4 * •— »» * ■ — ' 

Cj-Cj 2<?j * 2 2 J'll+l'l 


secondary leakage reflected 

load primary 


A *20 


T 2 

Note ti,»» K23+-J- just represents (with K3- Tj) the line to line admittance as expected. As before, 
the measred primary to secondary open circuit voltage transfer ratio 

~2.V[2 


8\i 


)'ir~yu 


A-22 


could be used to replaco the difficult to measure Similarly if the line to line secondary to open 

circuit primary voltage transfer ratio is where analytically 


£21 


”Tl2 

.Vll 


is meusured then 


t*| *« 


S 21 


1 *f 


A 

y\\ 


(C2-C3) 


A-23 


A-24 


and 


h 


„ , Yi y»i2 , ynKi\ 

__ r!J+T __ + .^_r 

j — — 

i'll 


A-25 


Other variations and measurements are possible. For example, if the following measurements were 
available 


y* 


Too “ 


h 


02 -Ci 


> , 22-y , 23 


(leakage effect) 


> 5 .— 00 

r ir r 2“° 


> r o 


y22-y:i _ yh 

2 >‘u 


C 2 ~Cj 

they could be used to obtain the desirable form 


A-26 


A-27 


Ij _ v„ , li , 

.. 1 2.1 + + y.u + w 

.Cj-Cj 2 . , >t 


A-28 


I + 


3'il 


NouaroumlciJ 
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Obviously lor this case, + 


is the linc-to-line secondary load and has been included to show that 


f V 

it does not effect the secondary to primary voltage transfer and enters into only as a parallel 

load, all as expected. The set of measurements y , h , r, , y lu seem appropriate for secondary to prt* 
mary in this case of floating line to line drive. Note that these measurements are redundant Uf V v 4 » is 
measured) and are related by 


Tot - Tu 't- Si\)'n 


A*28a 


u 
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Case A - Secondary to Primary: i\ m 0 % fj^O, Y$r* )%, )%— co 

This ease corresponds to a "line to neutral 4 * secondary drive with the other "line" free to take on a 
voltage determined by circuit considerations, It will he observed that Tor this case that the secondary to 
primary voltage transfer will be theoretically dependent on the connected secondary load ( \ ^ ) & 1 
Kv*). Note in order to show the dependence on the secondary loading Y\ has not been assumed equal 
to )%, Using equations A-2, A -4 (i|«?j*»0) obtain 

J p u0 , jj+^ , j2+ Yj) 


C i ** 


3®i2“*0’n4* ^i) 0‘jrh Yi+ Yi\) 


<?2 


3*12 

3'u4*i'i 


1 - 


y h +0'n+ Y\) i'jj) 


3'u~0'n4- )*i) )*2+ 1*23) 


l? 3 


A-28 


where the last form on the right is an algebraic rearrangement to show the difference from the line to 
line drive situation (see A-9 or A- IS). For reference 


C) « 


•1 + 


1*2 


yh 4- Oh rh F|) 0’3.\ — Yu) 
3'u~(3'ii 4- F|) 0^22+ fy4* ^ 3 ) 
h'h'~L v \& Y\) (3 ! 22“*3'234* 1*34*2 Fjy) ) 


<?3 


A -2 9 


A-30 


~~y h ”■* 0’ irl* l*i ) O’ 32+ )Vh Yu) 

where the last form on the right shows a difference factor from the "ideal" transformer situation/ 

% 

* , , . / 2 

The driving point admittance — can be obtained using the above and equation A-3, i,e. 

Ci 

h _ )'\i (2 ( 3*2 24*3* 23 ) 4- Fjj+bnjh Y\] b'lj — — 1*23(2 ( 3 * 224 * 3 * 23 ) 4 - Kj)— 

e l ~j2- y h — O'l 1 + l*t) 0*224* 1*3+ 1 * 23 ) 

secondary load shunting the ©j voltage drive 

It is of interest to observe: 

t'l 

a. The voltage transfer ~ (see A-28) depends on the secondary loading Y$ and Y 2i 

Cl 

b. If Yy-* 00 y j,e M shorting terminal 3 to ground yields (see A-28) the same voltage transfer ratio 

as A-18 or A- 19. No attempt so far has been made to examine the right hand form of A-2S 
c 1 

to examine how — varies with Y\ and Yjj, 
e 2 * ^ 

It is of interest also to note the special cases for Kj- F^«0, i.e. open secondary, i.e.. 




3 ' 1 2 0 s 234*3' 22 ) 

3* f2”" O' 1 1 4* F|)3'22 


<M 


3’ 12 


y n 4* Y\ 


1 


3®?24-0'||4* F|)3’ 23 


3* na — O' 1 1 4- F|)3 j 22 


Ci 


3'fi 4*0'||4* F|)3’23 
Cj — — ; — : rrr t’3 


1*2*" >V“ Yd^O 


3* f2 - 0’ u4* >' I )3'23 
h 0'234*.V22) [2y t f24*0'l|4* F|) 0'23”* 3*22)) 


Cy 


3'f2-0'l|4* 1*1)3*22 


A-31 

A-32 

A-33 


Because of the dependency on the secondary load which, of course, is highly variable, no attempt is 
made here to suggest a parameter model. 


* Nuio y parameters cannot be usoU to represent tin ideal U’.ti'isfonocr 
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Case B — Phase to Neutral Primary Connection, Neutral Not at Zero Potential 

When the neutral is not at zero potential, i.e., Y\ finite (see Figure 4), then certain modifications 
are required. In order to appreciate these modification requirements a short review oi' certain aspects of 
modeling is in order. First note that with a pole mounted distribution transformer it is obvious that the 
only current path to ground is via the intentional ground wire and associated ground rod, For vault 
mounted distribution transformers which obviously are in closer proximity of ground it is still reason- 
able to assme that any tank or conductor to ground stray capacity reactance is much greater than the 
intentional grounding impedance. The significance of these observations is that the whole transformer 
has a single known current path to ground (through a footing admittance Y.\) and the ellects of )\ are 
easily accounted for using the procedures discussed in Appendix l of Volume 3 “Technical Appen- 
dices” Final Report. Next note that for the same reasons the operation of the transformer (voltage 
transfers, driving point admittances, etc.) is unaffected essentially* by imperfect grounding. 

For outbound feeder to feeder model computed propagation, the only effect of the distribution 
transformer is to load the feeder from phase to neutral (for grounded wyS connection) or phase to 
phase .(delta connected system and for some cases of grounded wye systems). Letting Y t be the 
complex-scalar driving point admittance seen looking into the primary terminals of the distribution 
transformer the three phase admittance loadings are of the form 

0 0 0 
0 Y t 0 
0 0 0. 

Yr phase 2 connected 
^ >V-*~ 

for perfect grounding, i.e., F jV — oo, with imperfect grounding through a finite Kv these admittances 
become 


0 

0 

0 

0 


0 

0 

0 

0 

0 

Yt 

0 

-r, 


0 

Yt 

- Y t 

0 

0 

0 

0 

0 


0 

-Yj 

Y t 

0 

0 

-Y r 0 

Y t + r v 


0 

0 

0 

K.v, 


Y |- phase 2 connected A connected phase 2 + 3 

■=■ Ytfr^oa TjV^Oo 


The above are single examples of many combinations that can occur and serve to illustrate that given 
Yj which is unaffected by grounding considerations, the feeder loading admittance can be easily syn- 
thesized algorithmically for perfect or imperfect grounding. 


0 0 0 
0 Y t ~Y r 
0 - Y T Y t 

A connected phase 2 + 3 

y (V — eo 


• Conceivably a ground loop could exist vm the residential grounding. It is implicitly assumed that no voltage drop exists along 
the secondary neutral. 


r 


* 


; 4 

: ' i 
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3. Analysis of Simple KLC (plus ideal transformers) Two Winding Transformer 

This section examines a simple RLC (plus idea transformers) two winding model which takes into 
account struy capacity effects, as a single lumped capacitor across the primary terminals (a somewhat 
arbitrary choice). This is considered to bo a “zeroeth” order or very simplistic approach and is offered 
only to the extent that it may result in a model which can be used reasonably well to predict -.he 
response characteristics of two winding transformers at lower frequencies, say up to 50 kHz or possibly 
somewhat greater. A next higher order model might assume lumped capacitors between ail external 
terminals and from all external terminals to the case and/or core. Obviously as frequency range is 
further extended the ‘'finite element'’ representation would require further and further fineness to 
represent ali inter and intra-distributed capacity effects. 

The “zeroeth” order approximation is being suggested here for consideration because based on 
driving point admittance measurement made by R. Rankin 5/8/80 on a 25 kVA distribution 
transformer, it appeared that such a model could be used to fit the measurements quite well. However, 
no attempt was made to see how well it would predict voltage transfer ration and secondary to primary 
characteristics. 

Alsp it should be noted that this model is only applicable to: 

a. electrically symmetric distribution transformers 

b. balanced secondary loading 

c. L-L secondary drive 

d. phase to neutral primary connection 

The “zeroeth” order lumped parameter model is shown on Figure 2-1 and is observed to be the well- 
known 60 Hertz model with the lumped capacitor across the primary terminals. The above conditions 
(a, b, c) effectively allow the actual three-winding distribution transformer to be considered a two- 
winding transformer (since for these conditions no correct flows out/into the secondary mid-tap termi- 
nal). See “Magnetic Circuits and Transformers” (The Technology Press) for a good discussion of this 
model and especially leakage reactance. 
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Figure 2-1. “Zeroeth” Order Approximation for Two-Winding Transformer, with Stray Capacity Iron 
Core 


C 

Lx 

R\ 

Li 

R 2 

L m * 

R m ' 

N [t N 2 


a 


= “zerocth” approximation for stray capacity effects 
= primary leakage reactance effect 
= primary winding effective series resistance 
= secondary leakage reactance effect 
= secondary winding effective series resistance 

= equivalent series magnetizing inductance (can be frequency dependent) 
= equivalent series magnetizing loss (can be frequency dependent) 

= primary and secondary turns — ideal transformer 

N 2 


•The magnetizing shunt branch is usually, at lower frequencies, basically inductive with some dissipa- 
tion. Its representation could be either a series L m and R m as shown or a parallel shunt arrangement. 
In either representation the lumped parameters L and R will not be exactly frequency independent. 
For the series representation R is generally small and for the shunt representation R is usually large. 
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Tim two winding transformer will now be analyzed to derive various parameters identified in Part I. 
For convenience the analysts will be made in terms of notation of Figure 20 which has obvious rela- 
tionship to Figure 2*1. For Primary to Secondary we have; 

Open Circuit (Secondary) Primary to Secondary Voltage Transfer Ratio g^ 

l 


11 

Vi 


** Su *• 




a ( 1 + I’uZi) 

Open Circuit (Secondary) Primary Driving Point Admittance - 


II 


J--K.+ r ° 


z\\ H-VoZi 

Short Circuit (Secondary) Primary Driving Point Admittance y\\ 

_ y 

J ' U “ ' ;+ 2|+# J Z 1 (1+ )' U Z,) 

Primary to Secondary Voltage Transfer with Load 
t'i a 


v 2 


1 

» 1 

o (14* 1 oZ|) 

1 + 

ZrHrU* l\,Zi)Zi 



trU+loZi) 


Short Circuit (Secondary) Transadmittance y^ 


Jl 

vi 


- y n 


i'tr** 


Zi+n^ZjU+roZ,) 


Primary Driving Point Admittance — with Load 

Vi 


Jx y ll h + a+o 1 y»Zj) h 

n " L <i J (l+>\,Z,)+tZ r Wi i Z,U+>' l ,Z,)]>'i. 
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Figure 2-2. Two Winding Transformer Analysis Model 


Now this expression for the primary driving point admittance with load can be manipulated and 

v i 

expressed in terms of — , >*iii £121 T 12 as given above to vield 
-11 


ll 

‘’i 


•Vti + 



£|2 

^12 


Yl 
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Y 

which should be compared to equation A- 17a and where Yl is obviously equivalent to Kjj-f— ^ which 
is the line to line loading as assumed in Part A. 

Next we will look at the primary to secondary voltage transfer ratio. Equation 3-4 can be- manipu- 
lated to yield 


v± 

v\ 


1 

£12 
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which is exactly equivalent to equation A- 15 (realizing that the result in equation A- 15 needs to be 
multiplied by two to correspond to the line to line condition shown in 3-8). 

The purpose of the above exercise in deriving equations 3-7 and 3-8 has been to demonstrate that 
indeed the equations derived in Section 1 under the assumption of balanced secondary loading and 
symmetric distribution transformer electrical characteristics do represent a two winding transformer. 
We will not make any attempt now to derive the secondary (line to line) to primary characteristics. 

Reviewing frequency response data taken by J.J. LaForest it appears that the simple RLC-r ideal 
transformer being considered here* could account for most of the observed behavior out to 80-100 kHz, 
(.e., all the parameters seem to follow similar patterns. After 80-100 kHz a variety of “squirrelly” 
behavior is noticed due to probably the distributed-higher order stray capacity effects and possibly 
instrumentation vagaries. 

C. Conley and others have noted some additional resonance effects, usually small, in the 35-50 kHz 
region. Before adopting the simplified model suggested above, our transformer measurements should 
be designed and examined to see if other resonances occur. 

Assuming the simplified RLC + idcal transfoimcr model of Figure 2-1 is acceptable, then measure- 
ments can be designed to determine the appropriate lumped parameters. Generally it is not possible by 
measurements to determine separately L t and L ^ However, for most transformer loading situations, 
either on the secondary for primary tp secondary considerations or on the primary for secondary to pri- 
mary considerations, the current through the magnetizing admittance )» (Figure 2-2) is small compared 
to the current through Z], and therefore the IZ\ voltage drop is little affected by the current through 

To, At very light loadings, say open circuit, the series impedance effect ,Z|+-~* is dominated by 

H) • u 









i.e., Z\ 4* *rr ~ *4~ or H-F u Z,= l. 

* i) <11 

usually a reasonable approximation. 
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Net result of these considerations is that the following model is 



Figure 2-3 


Based on the preceding discussion concerning the observation that j *~~ | » | Z u !, i.e., 1+ Y Ki Z\ = 1,* 

Y o 

and examining equations 3-1 to 3-6 where the factor 1+ Y {i Z\ is seen to occur frequently, procedures 
for determining the lumped parameters C, L , n , R t , L\+a 1 Li > R\, R^ a can probably be developed l .\ith 
ingenuity. For example equation 3-1 indicates that the open circuit voltage transfer ratio should be rea- 
sonably constant with, frequency and equal to a~\ Equation 3-2 indicates, along with the model, that 

should asymtoticaily increase 6 db/octave at high frequencies and could be used to determine C 
2 u 

(see attached graph prepared by R.W. Rankin 5/8/80), Equation 3-5 or the model of figure indicates 
that y l2 could be used to determine R\+a 1 R 1 and R\ and R 2 can be frequently -estimated 

from design standards or approximated by dc measurements. Similarly % leakage impedance (primarily 
.reactance) data is available for all distribution transformers and can be correlated with measurements as 
appropriate. Other measurements such as y\\ (equation 3-3) could be used to help deduce the magnet- 
izing branch lumped parameters, either by a shunt or series representation whichever tends to lead to 
the least frequency sensitive components. Since L m and C tend to resonate at very low frequencies, 
approximately 300 Hz, such a resonant frequency test could be used to deduce L„ L m is usually more 
difficult to determine and various schemes (involving connecting the primary and secondary windings in 
series bucking and series riding are used and any book of transformer measurements could be con- 
sulted). However, at PLC frequencies, i.e., greater than say 3 or 4 kHz, we could ignore since it 
will be dominated by Y c , 

If the simplified model of Figure 2-3 is shown to be satisfactory, then an extremely simple model 
determination may be possible which will be good enough (since we never know exactly what the 
secondary loading is in practice, we probably cannot justify a high degree of perfection in the model 
anyway). The procedure would be to: 

1. Determine the effective capacity C using the asymtotic behavior of 

2 u 

2. Use nameplate available leakage impedance (reactance) specification to determine L^+a'Li 

3. Use known design Standards, to estimate R \ and Ri 

4. Ignore Y 0 

The adequacy of this approach should be given serious consideration providing our contractual upper 
frequency limit is reduced to 50 to 100 kHz. 

* For PLC rrequendes, i.e., 

V 7 R \+J"' L \ 
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APPENDIX 4 

PROGRAM FOR COMPUTING AND PLOTTING 
RLC TRANSFORMER MODEL PREDICTED RESPONSES 

R.C. Rustay and R.C. Wentz 


This Appendix describes a software program called ZPTRANS1 which given the RLC pa- 
rameters, computes, and plots the predicted bidirectional driving point admittance and voltage 
transfer ratio associated with the following transformer model. 


o 







This model is used to predict distribution transformer (DT) responses with the assump- 
tion that the Dt is electrically symmetric and has a balanced secondary load, (see Appendix 3) 
and where 

C — Lumped Stray Capacitance (Including Bushing Capacity) 

G * Core Conductance 
R -» Winding Resistance ^ _ . 

L - Winding (Leakage Inductance j Referred t0 P™ 1 "? s,de 
n — Voltage Ratio (essentially turns ratio for iron core) 

The model includes a frequency effect for the actual value of R. 

The attachment included with this Appendix is a verbatum copy of an internal memo writ- 
ten by R. Wentz further describing this program, its operation, and output. 
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DIFNAP SYSTEMS MEMO #52 
3/3/81 


Memo To: C, A. Stutt 

R. Wooding 
J. Gajjar 
R. Rankin 
R. Rustay 
R. Wentz 
A. Dunham 

From: R. C. Wentz 

Subject: New Program for Analysis of Two Winding Distribution 

Transformers 


Please be advised that there is a new program in the catalog that 
is available for use in the analysis of two winding distribution trans- 
formers. Output of the program includes a plot file which may be plotted 
on the Zeta plotter in Building #37. The source code of the program is 
contained in a file named ZPTRANS1 . 

The program utilizes a lumped parameter model to calculate four 
characteristics of the transformer, each over a range of frequencies 
from 2 kHz to 100 kHz. The four characteristics are: 

1) Voltage transfer ratio (in DB) 

2) Phase angle of the voltage transfer ratio 

(in degrees) 

3) Input admittance (in DB) 

4) Phase angle of the input admittance (in degrees) 

The plot file produced is composed of two separate graphs. The 
first is a plot of the voltage transfer ratio and its phase angle. The 
second is a plot of the input admittance and its phase angle. Also 
included in the output of ZPTRANS1 is an output file which contains the 
numerical values of the four plotted characteristics, the frequencies 
at which calculations were made, and the input parameters that produced 
those results. 

During operation of ZPTRANS1 , the user is requested to specify 
the following input parameters: 


1 ) Mode of Operation 

1 = Primary to Secondary 

2 * Secondary to Primary 

2) Load Characteristics - Used only for Mode 1 calcu- 
lations, The user specifies the values of load 
resistance and load inductance. For Mode 2 calcu- 
lations there is no consideration of loads "downstream" 
of the transformer. 

3) Feeder Characteristics - Used only for Mode 2 calcula- 
tions. The user specifies the values of a shunt feeder 
capacitance and a shunt feeder conductance. For Mode 1 
calculations, there is no consideration of feeders 
"upstream" of the transformer. 

4) Base Frequency - Since the winding resistance has been 
determined to be dependent on frequency, a base value 
must be known so that internal logic will produce the 
proper equation for frequency dependency. See me for 
details concerning this feature. 

5) Winding Characteristics - The user specifies values 
for the winding resistance and winding inductance, both 
referred to secondary. 

6) Turns ratio of the ideal transformer. 

7) Stray Capacitance - A shunt capacitance on the primary 
side of the model to account for stray capacitive effects. 

8) Core Characteristics - The user specifies values for the 
core inductance and core conductance. 

9) Output File - The user can specify the name for the 
output file or request output at the terminal, in which 
case no output file is written. 

10) Ordinate Range - The user can specify separate ranges 
for the ordinates of each of the four characteristics 
being plotted. 

Plots produced by ZPTRANS1 are intended to be on the exact same 
scale as those produced in the lab. This is to facilitate comparison 
of model predictions with results of laboratory measurements. 
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Attached is an example of a typical terminal execution of 
ZPTRANS1 , a sample page of the output file, and reduced copies of 
the two graphs produced by the Zeta plotter. 

If you have any questions concerning the operation of ZPTRANS1 , 
please do not hesitate to see me. 
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03/03/31 10.61 I 

THIS PROGRAM USES THE LUMPED PARAMETER MODEL 
TO ANALIZE DISTRIBUTION TRANSFORMERS 

MOCt 1 CALCULATES PRIMARY > SECONDARY 

MODE 2 CALCULATES SECONDARY — > PRIMARY 

ENTER 1 FOR MODE 1 OR 2 FDR MG D= 2 


PDF. cL-MENTS YOU mISH 10 IGNORE 


c.VTcR A VALUE OF 0. 


FEEDER CAPACITANCE = 
FtnDER CONDUCTANCE - 
cNThR DR IF OKAY, ELSE 


0. ' MICROFARADS 

0. 1 CC'GCC MHOS 
ENTER ! 


dASE F R EC u£NC Y FOR FREQUENCY DEPENDANT 
BINDING RESISTANCE = 20000. OCCCCC HERTZ 

cNTER CR IF OKAY, ELSE ENTER 1 


BINDING RESISTANCE = 
BINDING INDUCTANCE = • 
isNTER CR IF OKAY. ELSE 


0. I4C000 OHMS 
0.036 OCC A I LL1.HENRYS 
ENTER l 


* 


REFERRED TO 
SECONDARY 


TURNS RATIO Or IDEAL TRANSFORMER 

STRAY CAPACITANCE 

cNIER CR IF OKAY, ELSz ENTER I 


31.3 COOOO 

O.OCI 000 MICROFARADS 


EN1ER TURNS RATIO 


cNTER STRAY CAPACITANCE (IN MICROFARADS) 


4h$ 


Of 


CO 
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TURNS RATIO OF IDEAL TRANSFORMER 

STRAY CAPACITANCE 

ENTER CR IF OKAY, ELSc ENTER I 


CORE INDUCTANCE * 0. 

CORE CONDUCTANCE * 0. 

ENTER CR IF OKAY, ELSc ENTER I 


31.750000 

0. MICROFARADS 



MILLINENRYS ' 

MHOS 


* 

output file- *dtouti i 

ENTER CR rF OKAY, ELSE ENTER I 


DO YOU rtANT OUTPUT AT ThE TERMINAL ? 
cNTER CR rF OKAY, ELSE ENTER I 



ENTER 



E OF OUTPUT FILS 


OUTPUT FILE «DrOUT2 ? 

ENTER CR IF OKAY, ELSc ENTER I 

as 

J 

OPENED FILE * DT0Ur2 I 


i 


RANGE OF VTR MAGNITUDE * -30. 0 TO 
RANGE OF VTR FHASE ANGLES *-130.0 TO 
RANGc OF YIN-PUT MAGNITUDE * -50.0 TO 
RANGE’ OF YINKUr PHASE *-130.0 TO 
bNTER CR IF OKAY, ELSc. ENTER I 


20.0 D3 
160.0 DEGREES 
20.0 DB 
160.0 DeGRBSS 


* WROTE OUTPUT FILE * DT0UT2 * 

PLOTTING VTR’ VS. FREQUENCY 

PLOTTING VTR PHASE VS. FREQUENCY 

PLOTTING YINPUT VS. FREQUENCY 

PLOTTING YINPUT PHASE VS. FREQUENCY 

IN ROTE PLOT F I LE*RPLOTOO! ? 

i 

DO YOU MSH TO ANALIZE ANOTHER TRANSFORMER 1 
ENTER CR IF OKAY. ELSc ENTER I 


V 

* 


i 

is 
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. ********* ********************** ******$■* ****** ************************** 
OTOUTl 03/03/81 8J42 ANf 

****************************************** ******** *************** ****** 


hode of operation 
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BASE FREQUENCY FOR FREQUENCY DEPENDANT 

winding resistance = 20000.000000 hertz 


FEEDER CAPACITANCE 
FEEDER CONDUCTANCE 


0. microfarads 

0. 100000 MHOS 


winding resistance 

WINDING INDUCTANCE 


0.140000 OHMS REFERRED TO 

0.086000 MILLIHENRYS SECONDARY 


TURNS RATIO OF IDEAL TRANSFORMER 
STRAY CAPACITANCE 


31.750000 

0. 001000 MICROFARADS 


CORE INDUCTANCE 
CORE CONDUCTANCE 


0. MILLIHENRYS 

0 . mhos 


FREQ. (HZ) VTR (IN DB> 


VTR PHASE 


Y INPUT <DB) YINPUT PHASE 


cooo.o 

-1 .07 1 82 IE 

01 

-8 , 73724qE 

01 

-6.833358E- 

-01 

-8.736526E 

01 

2500.0 

-1 .265706E 

01 

-8.728262E 

01 

-2.6221H7E 

00 

*B.727362E 

01 

3000.0 

-1.424L45E 

01 

-8.716020E 

01 

-4.206578E 

00 

-8.716940E 

01 

3500.0 

-1 .55812 IE 

01 

-8.707449E 

01 

-5.546335E 

00 

-8.706189E 

01 

4000.0 

-1 .674I90E 

01 

-8.696929E 

01 

-6.707022E 

00 

-8.695489E 

01 

4500.0 

-1 .77658IE 

01 

-8.686627E 

01 

-7.730932E 

00 

-8.685007E 

01 

5000,0 

-1.868L82E 

01 

-8.676610E 

01 

-8.646947E 

00 

-8.674810E 

01 

5500.0 

-1 .951054E 

01 

-8.666902E 

01 

-9.475665E 

00 

-8.664922E 

01 

6000.0 

-2.026717E 

01 

-8.657506E 

01 

-1.023230E 

01 

-8.655346E 

01 

6500.0 

-2 .09632 7E 

01 

-8.648412E 

01 

-l .092840E 

01 

-8.646072E 

01 

7000.0 

-2. I60782E 

01 

-8.639608E 

01 

-1.157295E 

01 

-8.637088E 

01 

7500,0 

-2.220793E 

01 

-8.631078E 

01 

-1.2 17306E 

01 

-B.628378E 

01 

8000.0 

-2.276935E 

01 

-8.622806E 

01 

- 1 • 273447 E 

01 

-8.619926E 

01 

8500.0 

-2.329676E 

01 

-8.614776E 

01 

-1.326189E 

01 

-8.61 1716E 

01 

9000.0 

-2.379407E 

01 

-8.606974E 

01 

-1.375919E 

01 

-8.603733E 

01 

9500.0 

-2 .42645 IE 

01 

-6.599385E 

01 

-1.422964E 

01 

-8.595965E 

01 

10000.0 

-2.471086E 

01 

-8.591997E 

01 

-1.467598E 

01 

-8.588397E 

01 

10500.0 

-2.5 1 3546 E 

01 

-8.584799E 

01 

-1.510058E 

01 

-8.581019E 

01 

11000.0 

-2.554033E 

01 

-8.577778E 

01 

-1.550545E 

01 

-8.573818E 

01 

11500.0 

-2.592723E 

01 

-8.570925E 

01 

-1.589236E 

01 

-8.566785E 

01 

12000.0 

-2.629770E 

01 

-8 • 5642 3 IE 

01 

-1.6Z6282E 

01 

-8.55991 IE 

Oi 

12500.0 

-2.665306E 

01 

-8.557687E 

01 

-1.661819E 

01 

-8.553187E 

01 

13000.0 

-2. 699452E 

01 

-8.551285E 

01 

-1 .695964E 

01 

-8.546605E 

01 
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Figure (4-2) 
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APPENDIX 5 

MATRIX BASED GENERALIZED NEUTRAL REDUCTION PROGRAM 

R.C. Rustay 


Attached, as part of this Appendix, is a verbatum copy of an internal memo concerned 
with a generalized procedure to coalesce several electrically paralleled neutrals into one 
equivalent neutral so as to be compatible with a program FEEDPUSj. The requirement for 
this procedure arises usually from the analysis of multiphase cables as explained in the at- 
tached memo. 
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Appendix U 

Matrix Based Generalized Neutral Reduction Program 


The purpose of this Appendix is derive a general (matrix) based algorithm for the neutral 
reduction procedure required for the m*m, the Z and Y matrices encountered with under- 
ground cables. The Z and Y matrices are defined by the steady state coupled ordinary 
differential equations 


dV 

dX 

JL 

dX 


- Z / 

- Y V 


U-l 


U-2 


The neutral reduction procedure occurs because usually Z and Y are of a dimension,* 
say m, which is at least twice as large as the number of phases, n. This occurs because the ca- 
ble parameter programs consider also, as part of *he dimension m, the cable sheaths and any 
miscellaneous additional conductors which may be geometrically paralleling the cables, i.e., 
such as bare copper neutral wires. The purpose of the neutral reduction program is to obtain 
“reduced” matrices Z R and Y Ry each of dimension /i+l, by invoking the assumption that at 
every point position X , all of the /t-f 1 to m conductors (the conductor indexing is ordered so 
that the first n are associated with the phase conductors) are at equal potential K c , i.e., 


n+i 


Ki+2 “ 


= v„ 


U-3 


and that i e the equivalent current flowing through the equivalent (parallel) conductor is given 
by 


h * ffl+l "b hi+2 +••’•+* i m 
Thus given the representation of U-l and U-2, i.e., 


K 



Zj2 

v L 


Zl2 r 

^22. 

h 


‘ J'll 

^12 

h 


7i2 r 

^22. 


we seek reduced equations of the form 


K 


Zu 

CN 

N 

K 

am 

Z e 

z 22 


iu 


■^n 

Y\2e 

4 


. Y e 

yn , 


x 


x 


x 


x 


h 

h 

Vu 

Vl 



U-4 


U-5 


U-6 


U-7 


U-8 


where Z 12e and Y l2e are fl x l> Z e anc * Y e are an d Z 22 anc * ^22 are (complex) scalars. 


* Implied in this Appendix is the fact that Z, Y, Z R , Y R are square 
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Reduction of Y 

Equation U-2 can be partitioned to 


• 

h 


' I'll 

i'u 


k 

mm 

I’,/ 


X 


U-9 


or 


>V - I'll V U + I'll V L 

U-10 

>L - >V Y U + Yll Y L 

U-ll 

Introducing the m—//xl “summing" vector-matrix S 


[1 , 1 , • • • , 1] 

U-12 

then from U-4 


4-s r 4- s t y 22 v l 

U-l 3 

Now since all the components in V L are equal to see U-3, 


V L - 5 

U-l 4 

Thus U-10 and U-l 3 become 


i'u - i'll Y b - + K l2 S V, 

U-15 

4 - S r >V / (, + S T Y 21 Sv t 

U-16 


Equations U-15 and U-16 can be combined to form the desired reduced Kmatrix, i.e., 


>v 
h 

reduced Y matrix 


Ml 


r l2 s 

S T Y U S\ 


U-l 7 


It is dear that 


S T Y u T i s 1 

Y\ 2 S is m—nx 1 (transpose of above) 

S t Y 22 S is lxl scalar and simply sum of Y 22 elements 


Computationally the elements of Y n S could be formed by the individual row sums of J'n, 
and S t Y 22 S is simply the sum of all the elements, a scalar. Note that this algorithm holds 
for the case where «— 1, w«-2 (corresponding to a single phase cable). 

Reduction of 2 


Equation U-l can be partitioned to 




Z la 


X 


h 

4 


U-l 8 
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or 


i\-z n i L + z [2 i L 

U-19 

- Zr.% + Zsk 

U-20 

Inverting U-20 obtain 


4 -2a" 1 *■£ ~ z :: -1 W 

U-21 

Utilizing the "summing" vector-matrix of equation U-12 and proceeding as before (see U-13 
and U-14) obtain 




(. - 5 r Z,j-'S«V - S r Z 22 -'Z|2 r 4' 
Solving lor v e * noting S r Z., -1 S is n scalar, obtain 

. _ S^%£ _ i 

' s t z 22 -'s l s t z 22 -'s ‘ 

Next substitute U-21 into U-19 and obtain (using U-14) 

*V- ^\rZ n z 21 ^z n T) i v + z n z 22 -'s % , 


.Substituting U-23 into U-24 obtain 


n- 


Z 1 ,-Z,2Z J2 - | Z 1 , r + 


Z|;Z,;— l SS r 2;;- 1 Z,; r l 


s r z,,-‘s 


lr 4* 


Zi:Z22~ l S 

S T Z,~'S 


U-22 

U-23 


U-24 


U-25 


Now note that Z l2 Z 32 ~ l 5 is a single column vector-matrix and that S r Z 22 _ " 1 Z| 2 r «- 
(Zi 2 Z 2 2“' 1 5) r , since Z 22 ~ l is symmetric, and suggests a computational approach. Also it is 
observed that digital word precision could be a significant concern since an exact algebraic ap- 
proach would reveal many cancelling terms in the matrix coefficient of I u in U-25. 
Equation U-23 and U-25 can be combined to form the desired reduced Z matrix, i.e., 


K 


Z„- Z^Zvj^Zjj + 


z [2 z^ss r z^ z{ 2 
s T z~Js 


s r z^zj 2 

s T z~Js 



U-26 


reduced Z matrix 


Note that for the special case where w—1, m—2 (single cable situation), U-26 reduces to U-18 
as iit should. Also as before S T Z 22 ~ l S is simply the sum of the elements of Z 22 _1 . 

A rearrangement of equation U-25 using the identity 


Z 12 Z 22 “‘ l Zi2 


i Z} 2 Z22”*^5^'Z 2 2~^‘S'Z2 2 Z22~^ Z 
- S t Z 21 - [ S 


12 


U-27 


yields the alternate form 
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j 



Z 12 Z 2 - 2 ‘ [s%-jlsz 22 - ssr\ Zf 2 'Z 12 

1 

! z, 2 z 2 - 2 's 


z " s r z 2 - 2 's 

! s r z 2 - 2 's 

1 

1 


~r 

I 

S t Zu'Z { 2 

i 

i 1 

- 

S T Z$S 

1 STZx'S 

I 


where it is noted that SS T \s square matrix of dimension m-n with unity for every element. 
This upper left n*n submatrix in U-28 may be computationally more accurate than in U-26. 
(Remember 5 r Z 2 2 _1 5 r is a scalar.) 

Finally the comment is offered that the utilization of the m-nx 1 vector-matrix in the 
above analytical derivation does not imply that it needs to be used in the computational algo- 
rithm since by itself operating on a matrix represents, as appropriate summing the elements 
in each row or column. 
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